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ABSTRACT 



The potential energy of a static charge distribution on a lattice is rigorously computed in the stan- 
dard compact quantum electrodynamic model. The method used follows closely that of Weyl for ordinary 
quantum electrodynamics in continuous space-time. The potential energy of the static charge distribution 
is independent of temperature and can be calculated from the lattice version of Poisson's equation. It is the 
usual Coulomb potential. 
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The notion of approximating gauge field theories in continuous space-time by discrete lattice gauge mod- 
els grew (in part) with the desire to compute numerically the consequences of quantum gauge theories' 1-5 '. 
This might be accomplished without recourse to perturbative methods. However, some features of the dis- 
crete lattice version of QED theory appear more difficult than in the continuous space-time theory. Such 
appears to be the case for the force law between charges in compact lattice quantum electrodynamics' 6 " 13 '. 

For quantum electrodynamics in continuous space-time, the law of force between two static charges is 
that of Coulomb, while the law of force between two steady state current distributions follows from Ampere's 
law. For lattice compact quantum electrodynamics, the Ampere law can become somewhat difficult to 
compute due to the inherent nonlinear form of the field equations. On the other hand, the lattice compact 
quantum electrodynamic analog of the Coulomb law can be directly computed. This lattice Coulomb law 
calculation is the purpose of the work which follows. 

The compact quantum electrodynamic lattice model (to be discussed) is here described in the canonical 
Hamiltonian form. The cubic lattice of the theory (three space dimensions) has vertices located at positions 
{x}. The components of a vector field at a vertex live on the links {L} connecting neighboring vertices. In 
addition to the links connecting neighboring vertices, one must consider the plaquettes {P}. These are the 
square areas bounded by four links. An alternative view employs the dual lattice with vertices located at 
{x*} which are the center points of the cubes on the original lattice. The links on the dual lattice are in 
one to one correspondence to the plaquettes on the original lattice. Each link on the dual lattice punctures 
a unique plaquette area on the original lattice. 

The electric field components live on the links of the original lattice. The magetic field components live 
on the links of the dual lattice, or equivalently on the plaquettes of the original lattice as the Faraday law 
magnetic flux through the plaquette area. 

Perhaps the most simple manner to describe QED on a lattice is to employ a "quantum circuit array" 
viewpoint. (We use units where % = 1 and c = 1.) The (temporal gauge) relation between the electric field 
and the vector potential, 

E(x) = -A(x), (f) 
on the lattice appears as an expression for the voltage across a link 

V L = -4> L . (2). 

The energy stored in the electric field is defined by assigning (to each link) a capcitance C, so that the energy 
becomes 

K= {1/2)0^1 (3) 

L 

If we sum the voltages around the links (which form the boundary of a plaquette), then the lattice 
version of Faraday's law is that we must obtain the rate of change of the magnetic flux through the plaquette. 
Formally, Faraday's law on the lattice arises as a consequence of the definition of magetic flux $p through 
a plaquette. It is 

exp(iq$ P ) = Yl exp(iqt] L <j) L ) 7 (4) 

LedP 

where q is the "charge" unit of the theory, dP is the boundary of the plaquette, and the coefficients i]l = ±1 
dictate the orientaion of the plaqutte area. That we are discussing compact QED on the lattice now makes 
an appearance since the inductance A (assigned to each plaquette) represents a nonlinear circuit element; 
i.e. the magnetic field energy is given by 

U = (l/ q 2 A)Y,ll-cos(q<t> P )}. (5) 
p 

In total, the Lagrangian (L = K — U) of the model follows from Eqs.(2),(3) and (5) to be 

L = (l/2)Cj2ti (1/<Z 2 A) ]T[1 - coa(q* P )], (6) 

L P 
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together with Eq.(4). The implied Hamiltonian is 

H = -(1/2C) W> L ) 2 + (l/a 2 A) £[1 - «m(?$p)]. (7) 

L P 

The wave functions of the theory have coordinates that live on a circle and are thereby magnetic flux 
quantum (2u/q) periodic 

V(...,0L + (27r/ g ),...) = V(-,0L,...)- (8) 



The inner product 



< ijjf^i >= j ip* f (...,<f> L ,...)ipi(...,<f> L ,...)Y[{qd<l> L /2Tr), (9) 



completes the definition (quite rigorously for a finite lattice) of the Hilbcrt space of states. 

In the infinite lattice limit, there is one vertex per cube, three links per cube and three plaquettes per 
cube. On the other hand, there are only two linearly independent magnetic flux plaquette variables per cube. 
This arises because the magnetic field is transverse in the lattice sense to be described more fully below. 
Furthermore, if one starts at time zero in a subspace of states spanned by wave functions which depend only 
on the magetic flux coordinates 

= *(..., $ P ,...) (10) 

then one remains in this subspace at all future times. This is a consequence of the Hamiltonian Eqs.(4) and 
(7). This subspace corresponds to having no charges present on the vertices. The magnetic flux coordinates 
here completely determine the dynamics. Let us consider this in more detail. 
With link unit vectors defined as 

m = (1, 0, 0), n 2 = (0, 1, 0), n 3 = (0, 0, 1), (11a) 

each link on the lattice is determined by a vertex and a unit vector 

L^(x,n), V L =an--E(x) = (i/C)(d/d<f> L ), (lib) 

where a is the lattice spacing. The lattice divergence of the electric field at a vertex x is defined as 

3 

dwE(x) = (l/a)^n r [E(x) -E(x-an;)]. (12) 

i=i 

On the dual lattice we have the identities 

exp(iq$p) — exp[iqa 2 n ■ B(x*)], (13a) 

dw*B(x*) = 0. (136) 

The subspace of states spanned by functions of the magnetic field only, as in Eq.(10), are precisely those 
which obey the Gauss law with no charges at vertices 

dwE(x)^o = 0. (14) 

Now suppose that we place charges on a subset of vertices {x 3 }. The charge density (in terms of the 
discrete 5 function) is given by 

p(x) = (q/a^^z^x,^), (15a) 

3 

Zj = ±1,±2,±3,±4,... . (156) 
To compute the wave functions with the external charges we must solve the Gauss law equation 

dwE(x)V> = p(x)ip. (16) 
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To do so, we adopt the method used by Weyl^ 18 ) in continous space-time QED to the problem at hand, i.e., 
compact lattice QED^ 19 ). We first decompose the voltage accross a link into a transverse and longitudinal 
part; 

V L = -^(x,n)+x(x)-x(x + an), (17a) 

exp{iq(f)] = exp[iqcf>(x, n)] = exp[iaqn ■ A(x)], (176) 

divA(x) = 0, (17c) 

dwE(x) = {i/a 3 )[d/d X (x)}. (17d) 

Only the transverse link coordinates contribute to the Faraday law magnetic fux through a plaquette. Thus, 
Eq.(4) may be written in terms of transverse link coordinates as 

exp(iq$ P ) = Y\ exp(iq<f> L ). (18) 

LedP 

From Eqs.(17) and (6) we may write the Lagrangian as a "radiation" plus an "Coulomb" part 

L = (1/2)C^ t L - (l/2a 2 A) ^[1 - cos{0p)\ + L Coulomb , (19a) 

L P 
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Lcouiomb = (1/2)C ^ ^ [ X (x) - X (x + an, )] 2 . (196) 

X 2 — 1 

The Hamiltonian has a similar decomposition 

H ^radiation -^Coulombs (20(2 ) 

H radtatton = -(1/2C) ^(9/90 L ) 2 + (l/q 2 A) ^[1 - cos(q$ P )}, (206) 

L P 

Hcoulomb = -(1/2C) ]T G ( x - y)[d 2 /dx(x)d X (y)], (20c) 

where G(x — y) is the lattice Coulomb Green's function to be discussed in more detail below. Here we note 
that the subspace of wave functions of the form 

1> = exp[iqJ2 )]*(..., *p, ■..). ( 21 ) 

j 

solve the Gauss law Eq.(16) in virtue of Eqs.(10),(14),(15) and (17d). 

While the radiation problem in compact lattice QED is quite complicated, 

Hradiation^i-, $P, •••) = £*(•••, $P, .••). ( 22 ) 

radiation is decoupled from the Coulomb energy (for static charges if not for dynamical Fermions' 14-17 )). 
The Coulomb energy is simply 

(q 2 /2C)J2z i z j G(x i -x j ). (23) 

ij 

In going from the Lagrangian Eq.(19b) to the Hamiltonian Eq.(20c) one must compute the lattice Coulomb 
Greens function 

-V 2 G(x-y)=<5(x,y), (24a) 



in which the lattice Laplace operator is here defined as 
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V 2 /(x)-]T[2/(x)-/(x + an t )-/(: 



x - an,)]. 



(246) 



i=l 



Eq.(24a) may be solved by Fourier transformation as 



3 3 

G(x) = (1/2) I[f[(d9 t /2n)]e^^/[3-J2co S (e 3 )] 



(25) 




It is proved in the mathematical literature on lattice random walks that 



|x|G(x) - (a/47r) (|x| 



oo). 



(26) 



Our genral proof of the Coulomb law for compact lattice QED is thus completed. The physical motivation 
of the proof (which Weyl used in ordinary space-time QED) is simply as follows: If one enforces the Gauss 
law on the physical wave functions, then the Coulomb law will follow. That simple reasoning is sometimes 
obscured in the path integral formalism unless sufficient attention is paid to the lattice Helmholtz theorem. 

We make no statements as to nature of the free energy for the nonlinear radiation (transverse) degrees 
of freedom in lattice compact QED. But the Coulomb law for the longitudinal field components is on solid 
ground for all temperatures. In Fig.l, we compare the continuous space-time Coulomb potential to the 
Coulomb law obtained on an "infinite" lattice as in Eq.(25). The convergence to the continuous space-time 
Coulomb law takes about four lattice lengths. On the other hand, suppose that the calculation of G(x, y) 
were performed using a finite lattice with N x N x N lattice sites. At least N ~ 1, 000 would be required to 
obtain adequate numerical evidence^ 20 ) for a long range Coulomb potential. 
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FIGURE CAPTION 



Fig.l: The lattice Coulomb potential in the ni = (1,0,0) direction is compared with that obtained in the 
continuous space-time form of the Coulomb potential energy (q 2 Z\z<ij 'iirr). 
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